Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN KTpX style file v2.2) 



Absolute diffuse calibration of IRAC through mid-infrared and radio 
study of Hn regions 



^ '• Martin Cohen 1 *, Anne J. Green 2 , Marilyn R. Meade 3 , Brian Babler 3 

O ■ 

(N Remy Indebetouw 4 , Barbara A. Whitney 5 , Christer Watson 6 , Mark Wolfire 7 

-i— > ! 

q ! Mike J. Wolff 5 , John S. Mathis 3 , and Edward B. Churchwell 3 

ON 1 Radio Astronomy Laboratory, University of California, Berkeley, CA 94720, USA 
t-H | 

2 School of Physics, University of Sydney, NSW 2006, Australia 

r— I ■ 

^ 3 Dept. of Astronomy, University of Wisconsin, Madison, WI 53706, USA 

o\ : 

Astronomy Dept. University of Virginia, Charlottesville, VA 22904, USA 

m - 

O ' 5 Space Science Institute, Boulder, CO 80303, USA 

"Dept. of Physics, Manchester College, North Manchester, IN 46962, USA 

o ■ 

"^H 7 Dept. of Astronomy, University of Maryland, College Park, MD 20742, USA 

6 ■ 
a: 

^ i Accepted . Received ; in original form 



2 Martin Cohen et al. 



ABSTRACT 

We investigate the diffuse absolute calibration of the InfraRed Array Camera on the Spitzer 
Space Telescope at 8.0 /im using a sample of 43 Hll regions with a wide range of mor- 
phologies near /=312°. For each region we carefully measure sky-subtracted, point-source- 
subtracted, areally-integrated IRAC 8.0-/im fluxes and compare these with Midcourse Space 
experiment (MSX) 8.3-/im images at two different spatial resolutions, and with radio contin- 
uum maps. We determine an accurate median ratio of //MC8.0-/im/MSX8.3-/im fluxes, of 
1.55±0.15. From robust spectral energy distributions of these regions we conclude that the 
present 8.0-/im diffuse calibration of the SST is 36 percent too high compared with the MSX 
validated calibration, perhaps due to scattered light inside the camera. This is an independent 
confirmation of the result derived for the diffuse calibration of IRAC by the Spitzer Science 
Center (SSC). 

From regression analyses we find that 843-MHz radio fluxes of Hll regions and mid- 
infrared (MIR) fluxes are linearly related for MSX at 8.3 /im and Spitzer at 8.0 pm, confirm- 
ing the earlier MSX result by Cohen & Green. The median ratio of MIR/843-MHz diffuse 
continuum fluxes is 600 x smaller in nonthermal than thermal regions, making it a sharp 
discriminant. The ratios are largely independent of morphology up to a size of <~ 24'. We 
provide homogeneous radio and MIR morphologies for all sources. MIR morphology is not 
uniquely related to radio structure. Compact regions may have MIR filaments and/or diffuse 
haloes, perhaps infrared counterparts to weakly ionized radio haloes found around compact 
Hll regions. We offer two IRAC colour-colour plots as quantitative diagnostics of diffuse Hll 
regions. 

Key words: infrared: ISM - radio continuum: ISM - radiation mechanisms: thermal - radi- 
ation mechanisms: non-thermal - ISM: structure - ISM: supernova remnants 



1 INTRODUCTION 

Cohen & Green (2001: hereafter CG) made a detailed comparison 
of MIR and radio continuum imaging of Hll regions in an 8 deg 2 
Galactic field centred near I — 312°. Using 8.3-/u,m images from 
the survey of the Galactic Plane by the Midcourse Space experi- 
ment (MSX: Price et al. 2001) and 843-MHz continuum maps from 
the Molonglo Observatory Synthesis Telescope (MOST: Green et 
al. (1999)), CG developed a ratio of spatially integrated MIR to 
radio flux densities that discriminates between objects with ther- 
mal and nonthermal radio emission. This field was first explored 
at 843 MHz by Whiteoak, Cram & Large (1994: hereafter WCL) 
who characterized the emission of sources on all scales as thermal 



or nonthermal. The inventory of this region is diverse. WCL and 
CG list more than fifty thermal sources, ranging from compact Hll 
regions to filaments extending over 1° in length. Only about one 
third of the objects were detected in the recombination line survey 
of Caswell & Haynes (1987). Very few lie at tangent points so that 
no distance estimates, or only ambiguous ones, exist for most of 
the sources. 

The Spitzer Space Telescope (Werner et al. 2004: hereafter 
SST) GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Ex- 
traordinaire) program (Benjamin et al. 2003; Churchwell et al. 
2004) has surveyed 220 deg 2 of the plane in all four bands of the 
InfraRed Array Camera (IRAC: Fazio et al. 2004). GLIMPSE im- 
ages of the plane are roughly 1000 times deeper, and have an areal 



Diffuse 

resolution of order 100 times smaller, than those of MSX. The SST 
offers the prospect of exploring faint structure unseen by MSX. 

This study has several different objectives. All are achiev- 
able by combining radio and MIR continuum surveys of hundreds 
to thousands of square degrees of the Galactic Plane having the 
best resolutions currently available for such panoramic coverage. 
Spitzer provides GLIMPSE images of limited portions of the plane 
at high spatial resolution. First, we quantitatively investigate the 
diffuse calibration of IRAC using the photodissociation regions 
(PDRs) that envelop thermal radio regions. Our purpose is to at- 
tempt to link IRAC's 8.0-/im diffuse calibration to that achieved 
by MSX and validated absolutely to ~ 1 percent by Price et al. 
(2004). Absolute calibration was a vital activity on MSX, which 
had both point source and diffuse calibrations. Spatially integrat- 
ing images of point sources and comparing those results with the 
MSX point source catalogs shows these two calibrations are con- 
sistent. Therefore, MSX serves as a benchmark. Secondly, we con- 
firm CG's median ratio of MSX 8.3-^tm to 843-MHz flux densities, 
Fs.3 /5 , 843 , for thermal emission, using a broader assessment of sky 
background than CG used, and investigate whether this discrimi- 
nant varies with the density, angular scale, and morphology of ther- 
mal radio sources. Thirdly, we extend the MSX/radio discriminant 
between thermal and nonthermal radio emission to the equivalent 
IRACS.O-fim ratio, after allowance for the different radiance con- 
tributions of polycyclic aromatic hydrocarbon (PAH) emission in 
these two different filters. Fourthly, we examine IRAC images of 
compact Hll regions looking for evidence of MIR proxies for the 
weakly ionized radio haloes created by the leakage of photons from 
the dense cores. This phenomenon was proposed to explain the ex- 
istence of extended radio continuum emission associated with ul- 
tracompact Hll regions (e.g. Kurtz et al. 1999). Fifth, we compare 
the MIR spectral energy distributions (SEDs) of Hll regions in the 
field. Finally, we note that there are other radio structures in the 
field that have faint MIR counterparts that were not recognized by 
MSX because of its poorer resolution. One example is the opti- 
cally invisible planetary nebula PNG3 13.3+00.3, well-detected by 
GLIMPSE (Cohen et al. 2005). We draw attention to another par- 
ticularly interesting object and demonstrate that it too is a thermal 
radio emitter. 
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§2 introduces the theme of the absolute calibration of IRAC, 
explaining the origin of the differences between point source and 
diffuse calibration. In §3 we detail the 43 thermal objects stud- 
ied by CG (eliminating all sources known to be galaxies or to 
have nonthermal radio spectra). In §4 we describe our methods 
for assessing spatially integrated fluxes in MIR and radio images, 
and compare MSX8.3-/im integrated fluxes with the equivalent 
IRAC 8.0- fim values as a function of Hll region morphology. §5 
uses these SST/MSX ratios to probe the diffuse calibration of IRAC 
based on MIR spectra of our target Hll regions. §6 discusses the 
SEDs among the population of Hll regions in this field, expressed 
through their IRAC colours. §7 details an improved set of ratios 
of MSX 8.3-^m to MOST 843-MHz flux densities and provides 
the corresponding ratios constructed using the current calibration 
for GLIMPSE 8.0-/im image products. Statistical regression anal- 
yses are also given. §8 illustrates IRAC imagery of some of the 
regions studied, highlighting similarities and differences between 
MIR and radio morphologies. §9 presents a curious thermal object 
with unique characteristics. §10 provides a MIR/radio ratio for five 
galaxies and one supernova remnant (SNR) lying within this field, 
to illustrate the difficulty of measuring the MIR emission of SNRs 
even with the SST, but also to provide a quantitative estimate of 
the MIR/radio ratio for nonthermal emitters. § 1 1 offers our conclu- 
sions. 

2 THE ABSOLUTE CALIBRATION OF IRAC 

The basis of IRAC's absolute calibration is identical to the network 
of fiducial standards built by Cohen et al. (1999), founded on Sir- 
ius. This calibration has been validated absolutely by Price et al. 
(2004) to 1.0 percent. To support calibration with the high sensitiv- 
ity of IRAC required the development of new techniques and the 
establishment of faint optical-to-IR calibrators (Cohen et al. 2003). 
IRAC standards are either K0-M0III or A0-5V stars. The former are 
stellar supertemplates (UV-optical versions of the original cool gi- 
ant templates of Cohen et al. (1999)). The latter are represented by 
Kurucz photospheric spectra. These stars underpin the calibration 
of IRAC and provide cross-calibration pairwise between the three 
instruments of the SST. The primary IRAC suite consists of 1 1 stars 
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near the North Ecliptic Pole, but a secondary network close to the 
Ecliptic Plane offers efficient checks of calibration every 12 hr after 
IRAC data are downlinked. Reach et al. (2005) express their pref- 
erence for the A-stars and it is these that provide the fundamental 
parameters that are embedded in the header of each item of Ba- 
sic Calibrated Data (BCD) from the SST. Excellent and detailed 
summaries of IRAC's calibration are given by Fazio et al. (2004), 
Reach et al. (2005), and in the IRAC Data Handbook (version 3.0, 
available from the SSC.) 

Calibration of point sources is achieved by performing aper- 
ture photometry on BCD images, in a 10-pixel (12.2") radius aper- 
ture, and comparing with observations of the IRAC standards using 
the same process, aperture size, and sky annulus. Treating BCD 
images in this fashion results in consistent, well-calibrated, point 
source photometry that traces to MSX (Cohen et al. 2003). It had 
been assumed before the launch of the SST that the extended emis- 
sion calibration for IRAC could be derived from the point source 
calibration. Several tasks were planned to check the accuracy of 
the extended source calibration, but only after launch was it deter- 
mined that the responsivity to point source emission was lower than 
expected, and that differences existed between point source and 
extended source calibration, consistent with internal scattering of 
the light in the IRAC arrays. As a consequence, extended emission 
will not be correctly calibrated and must be corrected to the "infi- 
nite aperture" limit by multiplying by wavelength-dependent "ef- 
fective aperture correction" factors. The factor advocated for IRAC 
extended sources at 8.0 /im is 0.737 (IRAC Data Handbook, Ta- 
ble 5.7; Reach et al. 2005). The recommended factors for IRAC's 
three shorter wavelength bands are 0.944, 0.937, 0.772 at 3.6, 4.5, 
5.8 /im, respectively. In this paper we investigate the correction fac- 
tor at 8.0 nm by an independent approach. 

3 HII REGION INVENTORY OF THE FIELD 

Table 1 lists the thermal radio objects in the WCL study of the field. 
Thirteen sources listed by CG were excluded from our primary 
analysis, eleven of them nonthermal emission regions: five galax- 
ies, four known SNR, a further two nonthermal radio sources (one 
of them mixed inextricably with thermal emission), one source of 



unknown nature that is surely dominated by thermal dust emission 
(G313.28-0.33), and the source G31 1.31+0.60 (which is a mis- 
nomer by WCL for G3 1 1 .42+0.60 and appears in CG as a duplicate 
of this latter source's entry. G3 10.80— 0.38 cannot be isolated un- 
ambiguously in the complex Spitzer images). We have also inserted 
one new source, G3 11.70+0.31, of particular interest because of its 
unique morphology (§9). We initially adopted the source charac- 
ter (e.g. classical, compact, filamentary, etc.) from WCL. However, 
these attributes were drawn from heterogeneous literature so we de- 
cided to re-examine the Molonglo data for every source and have 
assigned our own homogeneous radio types to the sample. These 
are indicated in Table 1 (col.(3)), as are WCLs radio dimensions in 
arcmin (col. (2): Galactic longitude by latitude). Columns (5), (6), 
and (7) respectively present the observed ratios of the IRAC 8.0- 
fj,m to MSX8.3-/im flux densities, and of the MSX8.3-/im and 
IRAC 8.0- fim flux densities to those at 843 MHz. Uncertainties are 
given for each ratio (la) based upon the root-sum-square (RSS) of 
the fractional errors in the pairs of flux densities. Column (8) gives 
the areally integrated 8.0-fim SST fluxes in mjy. Column (9) gives 
figure numbers for objects represented by images. 

SST offers MIR images with far higher spatial resolution than 
those of currently available panoramic radio surveys. Therefore, we 
have also developed a MIR morphology based on the IRAC images 
(some of which are presented in §8), which appears in col.(4) of 
Table 1. Note that G3 11.64+0.97, a bubble some 35x30 arcmin, 
extends beyond the northern latitude coverage of GLIMPSE. For 
the calculations in this paper we have used the brightest portion 
that can be accommodated by the GLIMPSE latitude limit, namely 
the 6x 14 arcmin eastern filament centred at (31 1.88°, 1.01°). 

For convenience we have ordered sources in Table 1 by lon- 
gitude. To render the table more readable we have selected a small 
number of morphological descriptors that are coded by their lead- 
ing letters. In the MOST radio continuum these are: C - compact; S 
- shell (curved morphologies); B - blob (lacking any specific struc- 
ture but larger than a compact or unresolved source; and F - fila- 
mentary. The IRAC descriptors are: UC - ultracompact; C - com- 
pact; S - shell; D - diffuse; F - filamentary. However, it was neces- 
sary to add a second (lower-case) letter to signify secondary charac- 
teristics recognizable in the detailed imagery. These extra qualifiers 
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are: c - compact; m - multiple (when multiple compact cores occur 
within a source; s - shell; f - filaments present; h - halo; i - ob- 
jects that are apparently in close enough physical proximity to be 
interacting. 

The scheme we have adopted, based on GLIMPSE, must 
clearly represent those characteristics that are detected in this par- 
ticular survey. A deeper survey might well be expected to dis- 
cover fainter features unknown from GLIMPSE. A future change 
in pipeline, or in our understanding of diffuse calibration, could 
readily alter underlying background levels and local brightnesses, 
thereby affecting areal integrations of specific parts of an Hll re- 
gion. Similarly, longer exposures could reveal fainter structures 
that, at GLIMPSE'S depth, appear separate and distinct but, with 
greater depth in imaging, are found to be connected by filaments or 
to be merely small portions of widespread faint haloes. Therefore, 
we rejected the ratio of components' integrated fluxes as the deter- 
mining factor for morphological types. Instead we have assigned 
types by ranking the highest surface brightnesses found in continu- 
ous topological elements in residual images, typically at 8.0 /im. 

For example, Fig. 1201 (in §8) shows compact cores, a large, 
faint, partial shell, many separate, small, internal shells, filaments, 
and even diffuse emission. The highest surface brightness pixels 
occur in the very bright central core. The second highest levels are 
found in the small, complete rings north of the bright core, and 
the third in the long, curving filament that emerges to the north- 
west from the core, looping around to the north and almost reach- 
ing the brighter of the two small rings inside the large, faint oval 
shell. Hence, "Cs" is our designation. Tertiary characteristics are 
not given for simplicity so the secondary code describes whichever 
aspect is brightest after the primary. Had we used a third descriptor 
then G31 1.48+0.37 would have been "Csf. Parentheses around a 
primary code imply that it is not definite that the MIR source is the 
counterpart of the radio object. 

4 COMPARISON OF SST AND MSX INTEGRATED 
FLUXES 

CG used matched-resolution MSX and radio images to derive their 
MIR-to-radio ratios. The Galactic Plane is particularly rich in both 
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local and more distant structures which emit in the PAH bands. Ac- 
curate photometry of extended sources in the presence of this clut- 
ter requires special care. This affects MSX8.3-/im measurements 
and potentially three IRAC bands (3.6, 5.8, and 8.0 /im) containing 
PAH fluorescent features at 3.3, 6.2, 7.7, and 8.7 /im. The coarse- 
ness of the low-resolution MSX images could have caused us to dis- 
criminate poorly between faint MIR emission and sky background. 
For example, we might have sampled inappropriately bright regions 
to assess the sky emission, particularly if we now wish to seek faint 
outer haloes. Therefore, we have twice remeasured all regions in 
the images with 46" resolution (36" pixels), paying particular at- 
tention to the choice of local sky background and using more than 
a single local area to assess the variations in sky background when- 
ever possible. We have repeated this same procedure on the 20" 
resolution (6" pixel) "high-resolution" MSX images, better to as- 
sess the detectable extent of the sources, and with the same consid- 
erations of the sky emission to be subtracted. In this paper, we use 
the average of these two sets of measurements to provide our best 
estimates of the 8.3-/im flux densities of these extended regions. 
We have quantified the likely uncertainties in these integrals using 
the dispersions of the several sets of sky-subtracted measurements 
for each source. The average fractional uncertainty in these 8.3-^im 
areal integrations was found to be 6 percent. 

For the IRAC fluxes in all four bands we have likewise sought 
to include all IR emission apparently associated with the radio 
sources in our spatial integrals, and to identify background that is 
well outside the radio sources. Unlike the MSX images, we have 
worked from GLIMPSE "residual images". These are l.l°x0.8° 
images with 0.6" pixels from which all GLIMPSE point sources 
have been removed. The residual images are ideal for enhancing 
the recognition of diffuse nebulosity in regions of high point source 
density. They enable far more reliable photometry of such emission 
than can be obtained in the presence of the many point sources in 
the field. The residual images reflect the subtraction by our adapta- 
tion of DAOPHOT for GLIMPSE of all sources detected down to 
la, deeper than our publically released Catalog and Archive point 
source lists that extend down only to 10 and 5 sigma, respectively. 
Thus there are faint 2—3-a points that are subtracted from the resid- 
ual images that are not listed anywhere in our lists. Residual images 
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Table 1. The sample of Hll regions with radio and MIR morphologies, integrated MIR and radio flux ratios, and 8.0-^m fluxes 



Source 


Size 


Radio 


IRAC 


IRAC/ 


MSX/ 


IRAC/ 


IRAC 


Fig. 




(arcmin) 


morph. 


morph. 


MSX 


MOST 


MOST 


mJy 




G3 10.63-0.43 


<1 


C 


Cf 


1.64±0.20 


23±1 


37±5 


17600 




G3 10.68-0.44 


<1 


C 


Cf 


1.49±0.20 


23±1 


34±4 


13800 




G3 10.69-0.31 


<1 


c 


UC 


1.16±0.15 


5.0±0.2 


6±1 


5100 




G310.80-0.38 


7x6 


s 


Dc 




8±1 








G31O.89+0.01 


<1 


c 


Cs 


0.93±0.12 


37±2 


34±4 


7200 


18 


G3 10.90-0.38 


<1 


c 


Ch 


1.10±0.15 


11±1 


12±1 


9600 




G3 10.99+0.42 


5x7 


s 


S 


2.02±0.29 


120±10 


250±30 


363000 


23 


G311. 12-0.28 


7x8 


s 


Dc 


0.98±0.13 


15±1 


15±2 


50000 




G3 11.20+0.75 


6x15 


s 


Cf 


3.63±0.63 


46±6 


170±20 


325000 




G3 11.29-0.02 


15x14 


F 


Dc 


6.27±0.78 


40±4 


250±40 


831000 




G3 11.30+0.90 


8x4 


S 


D 


1.14±0.15 


72±2 


82±10 


64000 




G3 11.37+0.79 


8x24 


F 


Dc 


1.73±0.27 


87±9 


150±20 


298000 




G3 11.42+0.60 


<1 


C 


Cs 


1.83±0.23 


74±5 


140±20 


25000 


16 


G3 11.48+0.37 


5x8 


s 


Cs 


1.19±0.16 


24±1 


29±4 


129000 


20 


G3 11.50-0.48 


8x8 


s 


D 


2.16±0.29 


15±0.5 


32±4 


111000 


24 


G311.53-0.13 


2x9 


F 


Df 


1.93±0.26 


20±1 


38±5 


21000 




G31 1.59-0.59 


<1 


C 


Cf 


1.61±0.20 


17±1 


28±3 


7600 


17 


G3 11.62-0.60 


11x12 


s 


Dc 


1.31±0.18 


66±7 


86±10 


189000 




G3 11.63+0.29 


<1 


c 


Cm 


2.63±0.32 


13±1 


34±4 


32800 




G3 11.64+0.97 


35x30 a 


F 


F 


1.70±0.23 


48±5 


82±10 


43600 




G311.70+0.31 


3x2 


F 


C(§9) 


1.71±0.32 


27±5 


46±7 


4800 




G31 1.87-0.24 


<1 


C 


Ch 


1.61±0.21 


18±1 


29±3 


9000 




G3 11.90+0.08 


6x6 


S 


Cf 


2.30±0.31 


26±2 


59±7 


340000 




G311.93+0.21 


6x5 


s 


Cf 


1.64±0.22 


25±1 


42±5 


274000 




G3 11.97-0.05 


6x4 


B 


Dc 


2.30±0.46 


33±6 


75±11 


65700 




G312.1 1+0.31 


<1 


C 


Cs 


1.47±0.20 


28±1 


41±5 


42400 




G312.36-0.04 


<1 


C 


Df 


1.13±0.19 


24±3 


28±3 


13200 




G3 12.45+0.08 


60x90 


F 


F 


1.46±0.35 


800±200 


1200±230 


9460000 


26 


G3 12.50+0.32 


5x10 


F 


D 


1.56±0.24 


53±6 


83±10 


68500 




G3 12.58+0.40 


8x7 


B 


Cm 


1.85 ±0.24 


37±2 


69±9 


57000 




G3 12.59+0.22 


3x4 


B 


Df 


2.34±0.40 


29±4 


68±10 


45800 




LrJ IZ.OU+U.Uj 


< 1 




Lll 


o a Aj~r\ i 1 
z.4D±U.Jl 


AO 1 Q 


izU± i j 


13 /UU 


1 o 


G3 12.67-0. 12 


3x1 


B 


s 


1.46±0.24 


57±8 


83±12 


12300 


22 


G312.68+0.04 


3x4 


Ci 


Si 


2.39±0.35 


31±3 


75±10 


38300 


21 


G312.71+0.02 


3x3 


Ci 


Si 


1.05±0.15 


33±4 


35±5 


18000 


18 


G312.72-0.14 


<1 


C 


(S) 


0.98±0.13 


83±2 


82±10 


15900 




G3 12.77+0.06 


<1 


C 


Ch 


1.94±0.25 


250±20 


500±60 


28600 


15 


G312.93-0.09W 


8x9 


s 


D 


1.75±0.23 


59±4 


100±15 


88000 




G3 12.95-0.44 


5x4 


s 


Si 


1,06±0.18 


40±5 


43±5 


60900 




G313.07+0.32 


5x4 


s 


Df 


1.24±0.26 


64±15 


79±15 


30000 


25 


G3 13.32+0.05 


8x3 


F 


D 


0.98±0.13 


30±2 


21±3 


60600 




G3 13.37+0.02 


9x3 


s 


Dc 


1.27±0.17 


74±1 


94±11 


37700 




G313.46+0.19 


<1 


c 


Ch 


1.79±0.24 


16±1 


28±3 


22000 




Footnote: a We used the maximum bri 


ght area of 6x 14 arcmin of this region 


that was covered by GLIMPSE 



still contain sources not extracted by DAOPHOT, such as saturated 
sources and sources that peak beyond the non-linearity limit for 
each band. These objects were individually integrated and their flux 
densities subtracted from the integrals over the relevant regions. 
(Fig. I27l illustrates a residual image.) PSF-subtracted residual arte- 
facts consist of both negative and positive portions, often with a 
near-zero net integral. These usually have very little effect on in- 
tegrations of diffuse emission but they can make contributions of 
order 1 — 10 percent to small areas of integrated sky background. 
Therefore, we removed all such artefacts that lay within regions we 



took to assess the sky, and likewise subtracted the brightest recog- 
nizable artefacts lying within the boundaries of Hll regions. 

4.1 A statistical approach to removing MSX point sources 

To compare the truly diffuse component of emission for an Hll re- 
gion measured from MSX imagery with that assessed from the cor- 
responding IRAC residual image one must subtract the contribution 
made to the total MSX brightness by point sources. This is required 
to avoid biasing the direct ratio of IRAC and MSX brightnesses by 
the inclusion of faint point sources in the MSX images. We empha- 
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size that, were this the case, the ratio obtained would be lower than 
the true value because the MSX brightness would have been overes- 
timated. Ideally one should remove the identical point sources from 
MSX images as from the IRAC images but these are roughly 1000 x 
deeper than MSX so the vast bulk of GLIMPSE point sources were 
completely undetected, and undetectable, by MSX. 

The majority of GLIMPSE sources are normal stars that ap- 
pear faintest to IRAC in the 8.0-/im band, and were undetected by 
MSX. Any MSX point sources in the Hll regions that were dis- 
cernible to the eye (i.e. at least 3a) were integrated separately and 
their fluxes removed from that of the region. The median number 
of point sources subtracted per region was four. For a region with 
up to four point sources this approach was the simplest option and 
it was applied to more than half of our sample of Hll regions. For 
the seventeen regions having five or more point sources we queried 
the latest MSX Point Source Catalog (ver.2.3) using the "Gator" 
search tool at the Infrared Processing and Analysis Center (IPAC) 
and summed the flux densities of the resulting table of sources with 
detections within the specified area. 

In future investigations we would like to be able to remove 
the influences of point sources of MSX on wide area (tens to hun- 
dreds of square degrees) integrated fluxes as accurately as possible, 
including objects with signal-to-noise values of less than 3 — 5a. 
Such weak objects cannot be drawn from the MSX Reject Source 
Catalog due to the heterogeneity of this Catalog (that can compro- 
mise the validity of the sources themselves) and to the effects of 
Malmquist bias (that overestimate the flux densities for valid weak 
sources). Therefore, we have used the present paper as a formal 
test of a statistical method that relies on the use of the "SKY" 
model originally presented by Wainscoat et al. (1992) but with 
considerably expanded capabilities. Not only does the model re- 
liably predict the LogN-LogS distribution of bright and faint point 
sources in arbitrary directions, in any user-supplied system band 
from the far-ultraviolet to the MIR, it is also capable of summing 
up the smeared surface brightness due to unresolved and unrecog- 
nized point sources in any given area. This "Total Surface Bright- 
ness" mode was first described by Cohen (1993a), was applied to 
COBE/DIRBE diffuse sky maps (Cohen 1993b), and has under- 
gone testing since then in a variety of different spectral regimes 
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(e.g. Cohen 2000). Only one region in our sample contained so 
many MSX 8.3-um point sources that one could not remove them 
manually and one might not want to do so even by listing all sources 
in the MSX PSC2.3 and the Reject Source Catalog because of the 
likely impact of Malmquist bias on the resulting summation of flux. 
This area is G312.45+0.08, a 1.5° x 1° area containing the largest 
diffuse structures in this field and straddling the Galactic plane. It 
is, therefore, well-chosen in terms of sampling the highest surface 
densities of point sources as well as containing much clutter in the 
form of bright nebulous filaments. The newest version of the SKY 
model successfully predicts GLIMPSE IRAC source counts by in- 
cluding the warp of the disk and was used for the calculations be- 
low. 

The brightest and faintest stars in this field have magnitudes 
at 8.3 /im of —0.06 and 7.05, respectively (62 Jy and 88mJy). 
The statistical assessment of star counts in this 1.5 deg 2 area 
came from the computation that the total brightness of all point 
sources within these limits in the MSX 8.3-pim band amounts 
to 5.44±0.22x lO^ 11 W m- 2 deg~ 2 . Dividing by the MSX band- 
width and scaling by 1.5 for the total area observed yields 
582±23 Jy. We have tested these SKY predictions by directly find- 
ing all 8.3-pim point sources in the MSX PSC2.3 that lie within 
the rectangle defined by G3 12.45+0.08, by a corresponding J2000 
Equatorial polygon search using the IPAC Gator tool. Eliminating 
all sources that were undetected (i.e. have quality flag of or 1) 
results in a total of 693 point sources whose accumulated flux den- 
sity is 542±24 Jy, statistically the same (~ la) as SKY predicts. 
We conclude that SKY is a viable tool for assessing total surface 
brightness due to point sources in any region and we have utilized 
its capability in this paper to ensure that effectively the same set of 
point sources has been removed from the MSX images as from the 
IRAC residual images. 

SKY predicts more flux than that observed which is as ex- 
pected because the MSX Point Source Catalog may not be com- 
plete to 7.05 magnitude because of the clutter of the interstellar 
medium (ISM) due to PAH emission. Therefore, we also need a 
SKY calculation for MSX that mimics the DAOPHOT extraction 
process that generates the IRAC residual images, by removing all 
point sources to the 2a level, corresponding to 8.0 magnitude. That 
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calculation yields 625 ±26 Jy for the MSX point sources that would 
have been subtracted to the same brightness level as in our IRAC 
residual images at 8.0 fim. The observed sum of diffuse and point 
sources in the box integrated is 7100Jy. Therefore, point sources 
contribute only 9 percent of the integrated emission and the dif- 
fuse component is 6480 Jy in this very large area. For a more typi- 
cal diffuse region in our sample, with an area of ~90arcmin 2 , the 
point source contribution subtracted from the total integrated light 
at 8.3 /im is ~6.5 Jy, or 4 percent. Thirteen of the fifteen compact 
and ultracompact sources have three or fewer recognizable MSX 
point sources to magnitude 7.0, and their contribution is between 
0.3 and 2 percent of the dominating emission of the ionized zone 
itself. For the line-of-sight to the centre of each Hll region we 
ran the SKY model to compute the additional contribution made 
by point sources, fainter than the limit to which we could recog- 
nize sources in the MSX image, and extending down to magni- 
tude 8.0 at 8.3 to match the sources already removed from the 
coresponding IRAC residual image. The difference between assess- 
ing the MSX point source light extrapolated by SKY to magnitude 
8.0 and removing identified sources to magnitude 7.0 is only about 
7 percent of the contribution of point sources identified to magni- 
tude 7.0. Therefore, except for extremely large diffuse areas such as 
G3 12.45+0.08, in individual Hll regions the removal of recogniz- 
able MSX point sources accounts for about 1 — 4 percent of the total 
measured diffuse emission. The model extrapolations, to match the 
depth of removal of point sources in the IRAC residual maps, imply 
the subtraction of only about an extra 1 percent or less of the diffuse 
emission, well below the errors in the MIR integrated photometry 
of the same regions. 

We now summarize the corrections applied to the total inte- 
grated emission of regions for the contribution of point sources. 
The IRAC images from which we start the analysis are already 
point-source subtracted to a depth of 8.0 magnitude at 8.0 /im. We 
further remove artefacts by hand that correspond to slight mis- 
matches between our PSFs and real point sources. To subtract as 
closely as possible the same sources from MSX images as were 
removed in constructing the IRAC residual images ,we used a com- 
bination of techniques as follows, according to the character and 
apparent size of the Hll region. For compact and ultracompact ob- 



jects we measured fluxes of all recognizable point sources if fewer 
than five were seen in the image, and subtracted their total flux 
from the measured diffuse emission. For regions containing five or 
more point sources we extracted those using Gator and summed 
their fluxes. We extended every Hll region below the local 8.3-/xm 
magnitude limit for unambiguous identification in our MSX images 
by using the SKY-predicted total of surface brightness due to faint 
point sources in that direction. This contribution was subtracted 
from the photometry for each target region. We believe that this 
combination of techniques applied to the MSX images to account 
for point source removal is the best that can currently be achieved 
to match the depth of the IRAC residual images. If any starlight has 
been missed by the removal of actual sources and our extrapolation 
using SKY, it should be negligible compared with the diffuse flux 
attributable to any given Hll region. 

No statistical model of the sky, however physically realistic 
(as is SKY), can predict the occurrence of individual young star 
clusters either along the line-of-sight or within Hll regions. In- 
deed, objects in the bright cores of HII regions are observed on 
top of bright nebular emission and will not appear pointlike. Thus 
DAOPHOT will not identify them as point sources nor will they be 
removed from IRAC residual images. One might be concerned that 
such clusters would be treated differently by our handling of point 
source removal from MSX and IRAC images. However, our field 
contains no giant Hll regions nor complexes with multiple star- 
forming centers, commonly the sites of embedded clusters. None 
of our bright-cored regions, nor the more diffuse, open structured 
regions reveals any projected or apparent internal clusters to IRAC. 
This is largely because of the use of this field to demonstrate the ca- 
pabilities of the MOST survey by WCL. It is also a consequence of 
CG's choice to carry out their MIR-radio study, intended to focus 
on the ISM, and on interactions between the medium and individ- 
ually dominant stars, in this same field. Our analyses of MSX im- 
ages, therefore, suffer from no bias with respect to those of IRAC 
that is caused by unresolved star clusters. 

There are no low-resolution GLIMPSE images like the ones 
we used for MSX that offer quasi-independent photometry. Based 
on a series of tests carried out on diverse Hll regions of the repro- 
ducibility of our measured sky-subtracted fluxes we have assigned 
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Figure 3. Log-log plot of IRAC 8. 0-fim vs. MSX 8.3-/jm flux densities dis- 
tinguishing Hll regions of four different characters: ultracompact and com- 
pact; shells; diffuse; and filamentary. The key indicates the symbol associ- 
ated with the radio character of each object, based on our new MIR mor- 
phology. The two lines represent the regression relationships for 42 regions 
(dashed) and for the subsample of the 25 apparently smallest and densest 
sources (solid). 

a 12 percent uncertainty to the IRAC spatial flux integrals. The as- 
sociated tests involved varying the choice of local areas from which 
to assess the diffuse sky background. 

4.2 Direct comparison of MSX and IRAC images 

To illustrate the similarities and differences between MSX and 
IRAC images of the same regions we present grey-scale images 
of two fields (Figs. □ and |2}. MSX images of Hll regions are 
generally sparsely populated by point sources. We have chosen 
G3 10.99+0.42 because it shows a rather typical situation: many 
point sources but almost all of them undetected by MSX. The few 
that MSX saw are mostly avoidable by careful selection of the 
bounding rectangle for integration, or are readily isolated and in- 
tegrated separately (which we generally did). By contrast we also 
show a compact region apparently immersed in diffuse emission: 
G3 12.60+0.50, where SST shows a wealth of point sources but 
MSX does not. 

The fifth column in Table 1 gives the ratios of IRACS.0-/J,m 
to MSX8.3-/im integrated fluxes, with their statistical uncertain- 
ties. Fig.|3|presents the regression of these MIR fluxes, with errors 
in both variables. Symbols distinguish between Hll regions of four 
different characters as indicated in the key (ultracompact and com- 



pact; shells; diffuse; and filamentary.) The two lines differentiate 
between the regressions for the total sample of 42 regions avail- 
able with the SST (dashed) and for the subsample of the 25 regions 
with the best defined morphologies (solid); i.e. for all objects that 
we categorize as ultracompact, compact, or shells. The unusual ob- 
ject, G3 11.70+0.31 (§9), was included with these compact sources 
because of its small apparent size and well-defined structure. This 
subsample lacks all the diffuse and filamentary regions. 

Table 2 summarizes these results for each ratio of attributes 
by subgroup, giving the median, standard error of the median, size 
of the subsample, slope and offset with la errors for the regression 
between logarithmic quantities, and the Pearson correlation coeffi- 
cient between the two sets of MIR flux densities. The table is es- 
sentially divided between those regions that are well-defined and/or 
dense, those with complete or partial shell-like structures, and those 
consisting primarily of diffuse material, tenuous filaments, or bub- 
bles. Specifically the categories of Hll regions for which results 
are presented in Table 2 are: ultracompact and compact combined 
("U/C"); shells ("S"); ultracompact, compact and shells combined 
("U/C/S"); diffuse and filamentary combined (we have only two 
filamentary regions) ("D/F"); and the entire sample. Further conbi- 
nations of morphological subcategories reveals no meaningful dif- 
ferences in the ratio IRAC/MSX with respect to the character of the 
region. 

The regression for the entire sample has a slope of 1.07±0.03 
(1.00—1.14, ±3ct range), while that for the 25 densest and appar- 
ently smallest objects has a slope of 1.05±0.03 (0.98-1.13, ±3er 
range). There is no statistically significant difference between the 
two correlations. Both are consistent with linear proportionality be- 
tween SST and MSX 8-/im fluxes, have good correlation coeffi- 
cients, and span three and a half orders of magnitude in observed 
dynamic range of MIR flux. 

We have also carried out an independent analysis of the data 
offered in this paper using the types assigned by WCL (cometary, 
compact, classical, "Hn", "thermal", filamentary (including bub- 
bles), low-density). These are somewhat informal designations 
drawn from heterogeneous literature and there are no significant 
distinctions in terms of the ratio of IRAC to MSX 8-^im flux densi- 
ties even among WCL's types. 
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The entire sample available to define the ratio of IRAC/MSX 
near consists of 42 regions and yields a median of about 
1.6. There are no significant differences between this and any sub- 
sample of regions. The similarity of the ratios in physically very 
different sources such as diffuse and compact Hll regions con- 
firms the reality of this effect and marks it as a problem with the 
IRAC calibration. Our best quantitative estimate of the effect comes 
from an independent median analysis of the 25 smallest and dens- 
est regions: 1.55±0. 15. This is the empirical ratio of IRAC 8.0- 
Hm/MSX 8.3-/im found for these regions. 



5 MIR SPECTRA OF HII REGIONS: THE DIFFUSE 
CALIBRATION OF IRAC AT 8.0 MICRONS 

In addition to the SSC's formal study of aperture correction fac- 
tors, there is an analysis of the IRAC extended source calibration 
by Jarrett (2006), now an official SSC page 1 , based on a sam- 
ple of twelve elliptical galaxies whose observed light distributions 
have been modeled. That work finds a flux correction factor that 
varies from 0.78 to 0.74 at 8.0 fim, with a probable uncertainty of 
10 percent. This factor is the inverse of our IRACfMSX ratio and, 
therefore, implies ratios of 1.28—1.35, adopting MSX as the fidu- 
cial. That the ratios of IRAC/MSX exceed unity is believed to be 
due to scattering of light internal to the IRAC detector arrays, and 
it is a particular problem in the Si:As arrays (the 5.8 and 8.0-/im 
bands). The origin of the scattered light is thought to be in the 
epoxy layer between the detector and multiplexer. The glue has a 
strong spectral response, peaking in the 5— 6-/im region (Hora et 
al. 2004a). 

Not all of this excess is attributable to light scattering, how- 
ever. The IRAC 8.0-jj.m band has a very different relative re- 
sponse curve from the broader, somewhat longer, MSX 8.3-/im 
band (Fig.|4j, and the effects of this difference when integrating 
over the spectra of Hll regions must be treated. This we have done 
by three different methods. First we used the 2— 35-/im template 
for Hll regions, one of 87 such templates developed by Cohen 
(1993) to give the SKY model of the point source sky its capa- 

1 http://ssc.spitzer.caltech.edu/irac/calib/extcal/ 
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bility to predict source counts for arbitrary IR bands. This normal- 
ized template spectrum was constructed as the average of the ob- 
served spectra of 10 bright, rather compact (measured with beams 
between 3.4" and 30" in diameter) Hll regions. Within these IRAC 
and MSX bands the dominant spectral features from the cores 
(~1 arcmin in diameter) of Hll regions are (ranked by strength): a 
strong dust continuum, very deep 10-/^m silicate absorption; deep 
6.0-/im water absorption; PAH emission at 7.7 /im; and, some- 
times, weak fine structure lines (Willner et al. 1982). Integrating 
the two spaceborne bands over this (A,Fa) spectrum yields in-band 
fluxes of 1.75X10" 15 W cm -2 (SST) and 1.90xl0~ 15 W cm" 2 
(MSX). Dividing by the corresponding bandwidths (1.146 x 10 13 
and 1.402 x 10 13 Hz) gives 0.0153 and 0.0135 Jy, respectively. 
Consequently, in an optically perfect instrument, one would expect 
IRAC/MSX = 1.12±0.03. Away from the embedded central ioniz- 
ing source, the spectrum of extended emission from the PDR is still 
that of cool dust, overlaid by spatially extended silicate absorption, 
with diffuse PAH emission, and fine structure lines, even in large 
regions such as M17 (e.g. Kassis et al. 2002). Therefore, Cohen's 
template is still relevant because the bulk of the energy is emitted 
from the hot central core of an extended Hll region. 

The second method seeks to define the best average spectrum 
of the Hll regions in the field under study here, rather than relying 
on an average composed of the most famous and brightest compact 
regions. The Hll region template is a combination of ground-based 



Table 2. Analysis of ratios of /ft4C8.0-^im/MSX 8.3/im, MSX/MOST, and IRAC/MOST fluxes and sample sizes (in parentheses), and log-log regressions. 
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median± 


slope± 


offseti 
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coeff. 


sem(no.) 
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a 


coeff. 


sem(no.) 


a 


a 


coeff. 


U/C 


1.6±0.2(19) 


1.07±0.03 


-0.1 ±0.1 


0.97 


25±5(18) 


1.00±0.02 


1.3±0.1 


0.88 


36±13(18) 


0.98±0.03 


1.6±0.1 


0.82 


S 


1.5±0.3(6) 


1.14±0.07 


-0.5±0.3 


0.96 


45±20(6) 


1.11 ±0.05 


1.5±0.1 


0.87 


78±24(6) 


1.21±0.07 


1.3±0.2 


0.84 


u/c/s 


1.6±0.2(25) 


1.05±0.03 


-O.liO.l 


0.97 


27±13(24) 


0.99±0.01 


1.4±0.1 


0.83 


41± 14(24) 


0.93±0.02 


1.8±0.1 


0.79 


D/F 


1.6±0.4(17) 


1.11±0.04 


-0.3±0.2 


0.96 


44±18(18) 


0.77±0.03 


2.3±0.1 


0.72 


82±25(17) 


1.29±0.04 


0.9±0.1 


0.80 


All 


1.6±0.2(42) 


1.07±0.03 


-0.1 ±0.1 


0.97 


33± 18(43) 


0.93±0.02 


1.6±0.1 


0.77 


63 ±26(42) 


1.03±0.02 


1.7±0.1 


0.77 
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Figure 5. Calibrated LRS spectra for 1 1 Hll regions in the field. Error bars are ±lu at each wavelength. Our own, and WCL's (in parentheses), morphological 
types are given in each separate plot. 
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and airborne spectra, bridging gaps using data from the IRAS Low 
Resolution Spectrometer (LRS: Wildeman, Beintema & Wesselius 
1983). LRS coverage is from 7.7—22.7 /im and each spectrum con- 
sists of two sections, blue and red, whose effective apertures were 
6' x 5' and 6' x 7.5', respectively. For substantially extended 
sources the convolution of spatial structure and wavelength in these 
large apertures requires a complex extraction process. Sources 
smaller than 15" are not affected by the convolution (Assendorp et 
al. 1995). However, many Hll region spectra were published in the 
Atlas of LRS spectra (Olnon et al. 1986), based on a much simpler 
extraction method. This approach yielded perfectly adequate spec- 
tra (benchmarked by Cohen (1993) against ground-based and air- 
borne spectra in the construction of SKY's 87 templates) because 
most of the bright infrared emission of Hll regions arises in vol- 
umes of small apparent size. The complete Dutch LRS database 
comprised 171,000 LRS spectra, derived by using the same soft- 
ware as Olnon et al. (1986), and we have sought spectra for all the 
Hll regions in the / = 312° field in this database. We incorporate 
the recalibration of LRS spectral shape provided by Cohen et al. 
(1992a), and normalize the resulting spectra to the absolute 12-fim 
flux densities (scaled by the factor of 0.961 described by Cohen et 
al. (1996: their Table 6) as applicable to IRAS 12-/im flux densities 
of both the Point Source Catalog (PSC) and the Faint Source Sur- 
vey). When available we used the PSC flux density; otherwise that 
from the Faint Source Reject catalog. We found 93 LRS spectra for 
32 regions, and were able to extract useful averages for 1 1 of these 
regions. 

Fig. [5] presents the 1 1 averaged, calibrated spectra. Nine of 
these are somewhat similarly shaped; two others are quite distinct 
because of their blueness: G31 1.87—0.24 (Dc) and the cometary 
G3 12.77+0.06 (Ch). Each individual spectrum is plotted in Fa, and 
carries the source designation, and both our type and that assigned 
by WCL. Error bars are ±la at each wavelength. Quality varies 
from source to source although the characteristic deep 10-/im sil- 
icate absorption is seen in almost all the objects. Deep silicate ab- 
sorptions might indicate the presence of accretion disks that are 
viewed edge-on, but the frequency of these in Hll regions suggests 
that, in most regions, these absorptions occur because of deep em- 
bedding of the exciting stars in more spheroidal dusty envelopes. 



We created the average of the nine most similar spectra by 
inverse-variance weighting each spectrum (peak normalized to 
unity) based on the wavelength-dependent uncertainties that our 
reduction code produces (Cohen et al. 1992a). Fig. [6] compares 
this average LRS spectrum with the complete 2— 35-pim spectral 
template from the SKY model. The agreement between the two 
normalized spectra is very good between the LRS cut-on (7.7 fj,m) 
and the long wave limits of the bandpasses of the relevant MSX 
(e.g., 1 percent transmission at 1 1.08 /im) and SST filters (1 percent 
transmission at 9.63 /im). In the MIR, ionic lines do not contribute 
greatly to the content of the SST and MSX bands. 

We have constructed a hybrid spectrum by taking the data be- 
low 7.67 fim from the SKY template, and the average of the nine 
LRS spectra for our actual Hll regions beyond this wavelength. We 
integrated the MSX8.3-/im and SST 8.0-^tm relative spectral re- 
sponse curves over this spectrum, resulting in a ratio of the two 
instruments' flux densities of IRACMSX = 1.17±0.03. 

For the truly diffuse ISM the same calculations can be per- 
formed using spectra taken in the large-aperture (8x8 arcmin) 
near-infrared and MIR spectrometers of the Infrared Telescope in 
Space (IRTS: Onaka 2003, Sakon et al. 2004). We have assembled 
a spectrum from 1.4—11.67 based on the IRTS and for this the 
ratio of integrated flux densities, IRAC/MSX, is 1.15±0.06. The 
IRTS spectrometers were calibrated (Tanabe et al. 1997, Onaka et 
al. 2003, Murakami et al. 2004) with the identical basis to that em- 
bracing the IRAS LRS, MSX, and the SST so that integrations over 
this diffuse spectrum are directly comparable to those over the LRS 
and SKY template spectra. The same procedure can be applied to 
the ISOCAM CVF spectra of the ISM presented by Flagey et al. 
(2006). We computed the ratio IRAC/MSX for their spectrum with 
the most relevant longitude, the 299.7° field, and obtained 1.11. 
However, the traceability of the CAM CVF flux calibration is not 
identical to that of the MSX and SST and we assign this value 
twice the error of that derived from the IRTS spectrum. Conse- 
quently, whether we observe the most compact of cores, spatially 
extended MIR emission associated with Hll regions, or the ISM at 
large, we obtain a robust result for the measured ratio of IRAC/MSX 
near 8 /im with the empirical spectra of diffuse Hll regions, namely 
1.14±0.02 (by inverse- variance weighting of our four values). 
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Figure 6. Comparison between the SKY template for generic HII regions 
(solid line) and the average of 9 regions in this field (filled squares withilcr 
error bars). Open rectangles show the locations of the normalized median 
SED for the regions (§6). 
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6 SPECTRAL ENERGY DISTRIBUTIONS 



We have taken the spatially integrated flux densities (in Jy) from 
the images and multiplied by the bandwidths (in Hz) to gen- 
erate in-band fluxes. These we converted to Vega-based magni- 
tudes and, subsequently, into flux densities expressed in Fa (in 
W cm -2 /an -1 , after dividing by the bandwidths in /an. Table 3 
presents the absolute calibration for IRAC based upon the newest 
released relative spectral response curves (July 2004) and the ab- 
solute Vega spectrum used to define zero magnitude (Cohen et al. 
1992b) for IRAC and absolutely validated to 1.1 percent by Price 
et al. (2004). The assumed uncertainty on the response curves is 5 
percent, independent of wavelength. 



The actual ratio determined for IRAC/MSX is 1.55, of which 
the factor of 1. 14 is attributable solely to the different in-band fluxes 
measured by these two bands for the MIR spectrum of Hll re- 
gions. Internal scattering in IRAC's 8.0- /an detector, therefore, con- 
tributes the additional factor of 1.36±0.13 to the ratio of IRAC 8.0- 
/an/MSX 8.3-/tm based on our analysis of Hn regions. Expressed 
as the reciprocal quantity (the aperture or flux correction factor of 
SSC and Jarrett (2006)) our analysis is equivalent to 0.74±0.07. 
For comparison, the average correction factor and a conservative 
uncertainty of the elliptical galaxies at 8.0 /an is 0.74±0.07, while 
Reach et al. (2005) recommend 0.737 (no uncertainty specified) 
at 8.0 /an. Therefore, our independent analysis of the diffuse flux 
correction of IRAC at 8.0 /an based on Hn regions is consistent 
with two separate analyses based either on the performance of the 
instrument, or on models of the IR light distributions of elliptical 
galaxies. 

In the low surface brightness regime, the minimum detectable 
diffuse emission in this field has an intensity at 8.3 /an of about 
1.0xl0~ 6 W m~ 2 sr -1 . At the same location the 7ft4C8.0-/an 
level is 1.4 x 10~ 6 W m~ 2 sr -1 , consistent with the same factor for 
IRAC/MSX as measured for the bright Hll regions. 



It is necessary, however, to apply corrections for what is 
known about the shortcomings of the diffuse calibration of this in- 
strument. We have determined in this paper that the 8.0-/tm cal- 
ibration leads to fluxes too high by 36 percent. This can be com- 
pensated by applying an equivalent offset to the diffuse magnitude 
scale of +0.33. For IRAC's three other bands we adopt the SSC 
analysis (footnote in §5) based on model template SEDs for el- 
liptical galaxies. Averaging the "correction factors" derived from 
these galaxies indicates that corresponding offsets of +0.19, +0.09, 
+0.33 should be applied to measured diffuse magnitudes at 3.6 /an, 
4.5 /an, and 5.8 /an, respectively. Table 4 incorporates these correc- 
tions for each band both to magnitudes and to the values of Fa. The 
table summarizes integrated magnitudes of the sample of Hll re- 
gions giving the brightest and faintest magnitudes, and the medians 
with their standard errors. The dynamic range of the sample varies 
from a factor of 380 at 3.6 /tm to 200 at 8.0 /an, and the difference 
between the median magnitudes at these two wavelengths (almost 
5 mag) confirms the great redness of Hll regions. Every Hll region 
has its maximum IRAC F\ in the 8.0 /tm band. We normalized the 
quartet of Fa values for every object by dividing by Fa (8.0 /an), 
then took the median of each complete set of normalized values 
of FA(band)/FA(8.0/tm) where band = 3.6, 4.5, 5.8 /an. From the 
corrected magnitudes of the individual regions we have constructed 
Table 5, which presents the median values of the six colour indices 
between the four IRAC bands. Note that uncertainties are given for 
each offset based on the standard deviation of the average value 
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Table 3. IRAC zero magnitude absolute attributes 

Filter Bandwidth In-Band F^(iso) 

/im Won" 2 Wcm _2 /im _1 
Error Error Error 

(im % W cm -2 /im" 1 



with isophotal quantities and uncertainties below 

A(iso) Bandwdth F^(iso) i/(iso) 

fim Hz Jy Hz 

Error Error Error Error 

(im Hz Jy Hz 



IRAC1 


6.514E-01 


4.277E-15 


6.566E-15 


3.550E+00 


1.543E+13 


2.771E+02 


8.461E+13 




3.890E-03 


1.572E+00 


1.104E-16 


1.533E-02 


6.375E+10 


4.356E+00 


7.047E+11 


IRAC2 


8.827E-01 


2.342E-15 


2.654E-15 


4.494E+00 


1.306E+13 


1.794E+02 


6.685E+13 




5.707E-03 


1.592E+00 


4.560E-17 


1.994E-02 


5.834E+10 


2.853E+00 


5.657E+11 


IRAC3 


1.196E+00 


1.237E-15 


1.034E-15 


5.725E+00 


1.086E+13 


1.139E+02 


5.255E+13 




8.376E-03 


1.607E+00 


1.813E-17 


2.560E-02 


5.182E+10 


1.829E+00 


4.381E+11 


IRAC4 


2.395E+00 


7.230E-16 


3.019E-16 


7.837E+00 


1.146E+13 


6.310E+01 


3.865E+13 




1.663E-02 


1.599E+00 


5.263E-18 


3.528E-02 


5.242E+10 


1.006E+00 


3.172E+11 




2 4 6 8 10 12 

Wavelength (/um) 



Figure 7. Comparison of the median SED (±2 standard error rectangles) 
for thermal sources in this field with the diffuse ISM (solid line: IRTS, §4) 
and the SKY Hll region template spectrum (short dashed line). 

of the colour's correction factors for the nine galaxies. The errors 
in corrected colours are the RSS of the errors in the corresponding 
offsets and the standard errors of the medians. Those medians, with 
their standard errors, appear in col. (6) of Table 4. These normalized 
values represent a coarse SED over the IRAC range for Hll regions. 
Fig.Qdirectly compares this SED with the template for the cores 
of bright Hll regions and the diffuse ISM template. The SED is 
represented by rectangles, horizontally as wide as each bandwidth 
and vertically as tall as ±2 standard errors of the median. There is 
obviously much better agreement with diffuse emission than with 
the spectrum of a compact Hll region with its exciting star(s) and 
hot dust continuum. 

Figures|8|and|9|present two interpretive colour-colour planes 
for IRAC bands showing the 87 categories of source in the version 
of the SKY model where each source is represented by a complete 



2— 35-/im spectrum (Cohen 1993). The caption to Fig.|8|details the 
symbols used for various populations. These colour-colour plots 
can be used for quick-look diagnostics of populations in random 
Galactic directions. Planetaries are represented by two large filled 
circles because young and evolved planetary nebulae occupy two 
distinct regions in IRAS colour-colour planes (Walker & Cohen 
1988; Walker et al. 1989). Hora et al. (2004b) provide empirical 
IRAC colour-colour locations from their study of planetary nebu- 
lae with SST and their colours encompass the two larger circles 
predicted in Fig. [8] 

These figures show [3.6]-[4.5] vs. [5.8]-[8.0] and [4.5]-[5.8] 
vs. [3.6]-[8.0] diagrams. The large rectangles (±3 standard errors 
of the median) represent the zones in which diffuse Hll regions 
ought to fall based on the median corrected colours in col. (4) of Ta- 
ble 5. Correcting the observed colours using the offsets in col. (3) of 
Table 5 results in changes that are smaller than those that affect the 
magnitudes in any single band, and in no way alter the almost com- 
plete isolation of the diffuse Hll regions' box in these plots. Only 
one of the 87 types of point source in these plots encroaches on the 
diffuse colour box (the cross, corresponding to one type of reflec- 
tion nebula.) Fig. 1101 displays the observed colour-colour plots for 
the GLIMPSE Catalog of point sources for the roughly 1 deg 2 field 
centred at / = 312.50°. Ultracompact sources would pass into the 
GLIMPSE Point Source Catalog for this field provided that their 
cores were not saturated in any of the four bands nor contained 
any other corrupted pixels within the core radius. However, their 
colours would mimic the bright, compact Hll region colours repre- 
sented by the small filled circle in each of Figs.[8|and|9j It is inter- 
esting to note that there are three such pointlike sources observed 
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Table 4. Vega-based magnitudes of the sample of Hn regions in the IRAC bands, and normalized SED. 
Band X(pm) Brightest Faintest Medianis.e.med. Median F^(band)/F(8.0 fim) 



17 



3.6 
4.5 
5.8 
8.0 



2.24 
2.02 
-0.58 
-2.47 



8.68 
8.17 
5.21 
3.42 



5.60±0.25 
5.32±0.25 
2.71±0.25 
0.54±0.25 



0.23±0.05 
0.15±0.02 
0.66±0.02 
1.00±0.04 



Table 5. Observed and corrected median Vega-based colours (in magnitudes) of Hll regions in the four IRAC bands 



Colour 


Observed 


Offset±<7 


Corrected 




median 




medianis.e.med. 


[3.6] - [4.5] 


0.43 


+0.11 ±0.02 


0.53±0.06 


[3.6] - [5.8] 


3.35 


-0.13±0.03 


3.22±0.10 


[3.6] - [8.0] 


5.07 


-0.13±0.03 


4.94±0.10 


[4.5] - [5.8] 


2.91 


-0.24±0.03 


2.67±0.07 


[4.5] - [8.0] 


4.69 


-0.28±0.03 


4.41±0.08 


[5.8] - [8.0] 


1.79 


+0.00±0.04 


1.79±0.05 



with the colours of (ultra)compact Hll regions. There is only triv- 
ial contamination, by any of the 213563 observed point sources in 
the GLIMPSE Catalog for the 1 deg 2 region centred on 312.50°, of 
the corresponding colour-colour box that we have determined for 
diffuse Hll regions. 

The overall dispersion in SEDs expressed by the IRAC colours 
is quite small, in part because of the rather limited contributions of 
emission lines in most Hll regions. The primary features in their 
MIR spectra are attributable to broad absorptions by silicates, ices, 
and organic materials. Cohen et al. (2005) considered the emission 
in the IRAC bands from ionized gas in a new planetary nebula found 
in GLIMPSE observations. Both free-free and hydrogen line emis- 
sion were derived from the level of thermal radio emission but did 
not amount to more than a few percent in the IRAC bandpasses. 
Flagey et al. (2006) perform the same kind of estimates for the 
diffuse ISM, deriving no more than about a 12 percent contribu- 
tion from recombination lines and free-free continuum even at the 
shorter IRAC wavelengths. 

Fine structure lines are not bright in Hll regions as evidenced 
by their absence even in well-known compact sources (e.g. Win- 
ner et al. 1982). One can similarly gauge the minimal influence of 
emission lines in the LRS range through Fig. [5] Even for the IR- 
brighter sources in our field one sees only three 12.8-/^m [Nell] 
lines and one object might show an 18.7-/im [SIII] line. 




0.0 0.5 1.0 1.5 2.0 
IRAC3N-IRAC4N 

Figure 8. Diagnostic IRAC colour-colour plot ([3.6]-[4.5] vs. [5.8]-[8.0]) 
comprising 87 types of IR point source. Key: pluses - normal dwarfs, gi- 
ants, supergiants; asterisks - AGB M stars; diamonds - AGB visible C stars; 
triangles - AGB deeply embedded IR C stars; squares - hyperluminous ob- 
jects (these objects include deeply embedded OH/IR stars and early-type 
hypergiants (Cohen 1993; a small number are required to reproduce MIR 
source counts at low latitude: Wainscoat et al. (1992)); crosses - exotica (T 
Tau stars, reflection nebulae); larger filled circles - planetary nebulae; small 
filled circle - bright compact Hll regions. The reddening vector correspond- 
ing to an Ay of 10 mag is shown by the shaft of the arrow in the upper left 
corner. The large rectangle denotes the zone occupied by the extended Hll 
regions in this paper (±3 standard errors). Axis labels include magnitudes 
with "N" in their names to signify that new relative spectral response files 
from July 2004 were used for the IRAC bands. 
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Figure 9. As Fig.^but for [4.5]-[5.8] vs. [3.6]-[8.0]. 
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-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
[5.8]-[8.0] 

Figure 10. Same axes as in Fig. |9] but based on the 213563 observed 
GLIMPSE point sources in the catalogue for this field. The densest regions 
of the plot are represented by contours rather than showing every source. 
The appropriate Hll region box shown in Fig.|8]also appears in this plot. 



7 COMPARISON OF MIR AND RADIO FLUX 
DENSITIES 

Fie. II ll compares the MSX8.3-/im flux densities with those mea- 
sured from MOST images (both in mjy). Fig. ^| similarly com- 
pares IRAC&.O-fim flux densities with 843-MHz data. These two 
figures retain the different symbols to identify different types of 
thermal region. Table 2 summarizes the median analyses and re- 



gression lines for the subsamples of sources in terms of the ratios 
of MSX to MOST (centre of the table), and SST to MOST (right 
section of the table), flux densities. 

There are meaningful correlations implying linear proportion- 
ality between both MSX and SST 8-^im fluxes and radio contin- 
uum fluxes. However, unlike the IRACfMSX correlation, there are 
three issues which affect the MIR to radio flux ratios. One expects 
the MOST to lose flux substantially as the interferometer observes 
larger regions. The best model for the actual synthesized beam of 
the MOST incorporates an effective minimum spacing of about 70A 
(Hogan 1992) although this parameter varies between observations 
and may be dependent on meridian distance. Therefore, angular 
scales greater than about 30 arcmin are typically poorly detected. 
This certainly affects the largest structure in the field (the filamen- 
tary object, G3 12.45+0.08, MSX/MOST of 800) with a size greater 
than 1°. It might also overestimate MSX/MOST for the second 
largest source (the diffuse G3 1 1.37+0.79, MSX/MOST of 87), with 
a size of 24 arcmin, corresponding to the second largest Hn region 
in our whole sample. To assess the fraction of true radio continuum 
flux that has been measured for each specific region by the MOST 
would require detailed modeling based on the irregular structures 
of real sources, and on the u-v coverage obtained, like the simula- 
tions performed by Bock & Gaensler (2005). We can state only that 
MSX/MOST and IRACMOST ratios are invalid for G3 12.45+0.08 
but we suggest that this phenomenon may also contribute to some 
of the other large ratios that are found in Table 2, particularly for 
diffuse regions. 

The second concern is relevant to the single example of a 
cometary Hn region in our sample, G3 12.77+0.06. This appears 
anomalous in both Figs. 1 1 1 l and l 12l because it is quite compact yet 
also has a very large MIR/radio ratio. We believe this object is op- 
tically thick at 843 MHz, greatly diminishing the radio continuum 
flux density. This would be consistent with the extremely deep 10- 
fim silicate absorption whose optical depth, r(10/im) ~3, hence 
Av^55 (Roche & Aitken 1984; Chiar & Tielens 2006), and the 
fact that this source has the most extreme LRS energy distribution 
in our sample in that it shows the deepest 18-/im silicate absorp- 
tion too. The blueness of the spectrum implies an abundance of 
warm (>400 K) dust (rather than cool, i.e. <100 K) local to the 
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embedded high-mass star. Consequently, this region's exciting star 
is enormously embedded in dust, plausibly leading to local free- 
free absorption from the associated dense ionized gas. Therefore, 
we have chosen to exclude this region from any MIR/radio analysis 
of subsamples, reducing the ultracompact and compact Hll regions 
to 18 objects. However, samples designated as "AH" in Table 2 in- 
clude this cometary region. 

Thirdly, our SST images have been culled of point sources 
(only their residuals after PSF subtraction remain) but we have no 
such products for MSX. Thus MSX fluxes within the radio source 
rectangles could potentially still include some faint point sources 
although efforts were made to remove their influence. Radio flux 
densities for big sources will suffer loss of signal because of resolv- 
ing out structure as discussed above. Spitzer images are so sensitive 
to diffuse emission at low surface brightness levels that it could be 
necessary to seek uncontaminated sky at greater distances from the 
source than for MSX. This could lead to underestimated sky back- 
grounds, hence somewhat enlarged source fluxes for IRAC, and rel- 
atively reduced fluxes for MSX. Both effects could inflate the ratio 
of IRAC/MOST. Consequently, we chose matching sky regions for 
both IRAC and MSX to mitigate this problem. Our best estimate 
for these ratios will clearly come from the sample of the 18 most 
compact entities, excluding the cometary source. From this set we 
find MSX/MOST has a median of 25±5, in excellent accord with 
CG's median of 24. The corresponding empirical ratio (i.e. with no 
flux correction factor applied) for SST is 36±13. 

Figures fT3l and fT4l seek correlations between MIR/radio flux 
densities and the type of Hll region, and its characteristic size, 
respectively. These figures are derived from the data presented 
in Table 2. It is sufficient to explore any trend with the ratio 
MSX/MOST, because we know SST and MSX fluxes are linearly 
correlated. A weak trend toward larger median ratios in less com- 
pact sources is suggested by Fig. 1131 and explored in Fig. 1141 in 
which we have attempted to quantify the size of the different mor- 
phology types, using the square root of the median of the areas 
in each category. For ultracompact and compact regions we have 
measured the full width at half maximum (above the surrounding 
sky) of the bright cores at 8.0 and have used this as the char- 
acteristic dimension for these objects. The median size of these 18 
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Figure 11. Log-log plot of MSX 8.3-jim vs. S(843 MHz) flux densities. 
Symbols and regression lines as in Fig. [3] The compact object near the left 
boundary in the plot whose filled symbol is also encircled represents the 
sole cometary region, G3 12.77+0.06, which has been excluded from the 
regression analyses of compact regions. 
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Figure 12. Log-log plot of IRAC$.Q-fim vs. S(843 MHz) flux densities. 
Symbols and regression lines as in Fig. [3] 



compact regions is 25±4arcsec (standard error of the median), or 
0.42arcmin. The 15 regions whose ratio of MSX/MOST exceeds 
the median plus 3 standard errors of the median (a value of 40) 
consist of 3 compact, 2 shell, 7 diffuse, and 2 filamentary sources 
(excluding the self-absorbed G3 12.77+0.06). Each of these 3 com- 
pact regions also contains filaments, or shell structure or has a MIR 
halo. Therefore, this trend probably reflects the onset of the loss of 
radio signal with increasing apparent size. 

It might instead indicate that compact regions other than 
G3 12.77+0.06 are dense enough to suffer self-absorption of free- 
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free emission. Alternatively the trend might arise because the high 
density of dust in physically small sources leads to stronger MIR 
thermal continuum emission at the expense of fluorescent PAH 
bands. However, this would also increase the dust emission from 
small grains (PAH clusters with hundreds of carbon atoms) ob- 
served as a broad plateau of emission in the 8.0-/im IRAC band, 
underlying the 6.2, 7.7, and 8.7-/im bands. This would offset some 
of the lost PAH emission. Detailed modeling of specific regions is 
required to explore this alternative hypothesis, incorporating both 
the formation and destruction rates of PAHs and small grains, but 
this is beyond the scope of this paper. 

There is significant scatter in Fies. II ll and ll2l but one might 
have predicted even more than is observed because variations in 
excitation parameter and in 843-MHz optical depth from source 
to source would affect the radio spectral shape and turnover fre- 
quency. However, as CG concluded, most of the Hll regions in 
this field are not optically thick at 843 MHz; even the compact re- 
gions. The very fact that MOST detected so many regions argues 
that the old criterion, that below 1 GHz Hll regions are optically 
thick, is not sufficient. The general absence of optically thick re- 
gions at 843 MHz in the 1=312° field is also borne out by Figs. 1131 
and !14l The small MIR/radio ratios among compact objects implies 
there is little loss of 843-MHz signal. By comparison with the me- 
dian of 25 for MSX/MOST for compact regions plotted in these 
figures, the cometary G3 12.77+0.06 has a ratio of 250, illustrating 
plausibly how these ratios would be affected by significant optical 
depth. 

Among the many Hll regions in this field is one object that 
would be defined as ultracompact, namely G3 10.69— 0.31. The 
SST 8.0-/im image has all the characteristics of a minimally re- 
solved bright source; a size between 1.5" and 2" (slightly in excess 
of the point spread function's FWHM) and the hexagonal pattern of 
SST diffraction spikes. The residual image (with the suitably scaled 
PSF subtracted) shows a central hole at the core but has a saturated 
ring immediately around this hole to a radius of roughly 1 " . Its ratio 
of IRACfMSX is 1.16±0.08, exactly as predicted for a source with 
the SED of an Hll region that we have presented in Fig.0 This re- 
sult also confirms that IRAC photometry of a region dominated by a 
pointlike core (although we cannot exclude the existence of a very 
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Figure 14. As Fig. ll3l but the ratios are plotted against the logarithm of the 
median apparent size of objects in each morphological group. 



faint halo of emission) is essentially correctly assessed. Inspection 
of the Basic Calibrated Data images sugests that the object is at or 
near saturation at 8.0 /im, but not severely saturated. Therefore, we 
believe there is no major loss of flux in integrating its image and 
the IRAC/MSX ratio is credible. However, the source is too faint 
to have been detected by the LRS instrument. The ratios of both 
MSX/MOST (5) and IRACMOST (6) are the lowest of any source 
in Table 1, perhaps because dust competes more effectively for stel- 
lar UV photons than do PAHs. 
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8 MIR IMAGERY OF HII REGIONS 

We present monochrome SST images of twelve objects that span 
the range of MIR morphologies of Hll regions in this field. These 
are chosen to highlight the structural MIR characteristics that can 
be associated with various radio types. Although the printed paper 
presents black and white renditions of the IRAC images, the elec- 
tronic version shows them as 3-colour pictures in which 4.5-fim 
(blue), 5.8-/im (green), and S.O-fim (red) images have been com- 
bined. Our reasons for this particular trio are as follows. Firstly, 
3.6 /im shows too many stars and clutters the diffuse structure. 
MSX was often unable to provide the requisite confirmation of 
PAHs by the test of matching morphologies in its 8.3 and 12.1- 
um bands (both of which include PAH band emission) because the 
sensitivities of all its other bands were so much lower than that of 
the 8.3-/im band. By contrast, IRAC offers two sensitive bands ca- 
pable of sampling PAH emission. In order to confirm the existence 
of PAHs, that are widely believed to be the dominant spectral fea- 
tures in diffuse Hll regions, we wished to see the same morphology 
emerge in IRAC's 5.8 and 8.0-^im bands. One can also seek subtle 
changes in PAH populations through spatial variations in intensity 
ratios or, equivalently, in false colour across regions. For example, 
PAH emission from PDRs generally appears yellow or orange in 
our false-colour gallery because both 5.8 and 8.0-/im bands contain 
comparably strong fluorescent PAH emission features. With IRAC 
we still have a probe for H2 lines through the 4.5-/im band too. 

Although colour differences across a single Hll region have 
relative significance, note that variations between colour images 
of different regions do not necessarily indicate any spectral differ- 
ences between the objects. Our intent was simply to offer colour 
images that would reveal the intrinsic structure of each region, set- 
ting the colour balance and contrast object by object, rather than 
aiming at an absolutely black "sky" for all sources, or a yellow 
level for PAHs that is constant for every source. 

The MIR appearance of the cometary region in our sample 
(G3 12.77+0.06, Ch: Fig. 1151 clearly resembles the canonical radio 
structure yet it is also associated with a considerable number of 
filaments. 

The term "compact Hll region" may be something of a mis- 



- 




G312.77+0.06 



Figure 15. Cometary Hll region G3 12.77+0.06 (MIR morphology type Ch). 
with a diffuse halo. The image is orientated in Galactic coordinates and 
measures 205 by 228 arcsec, I by b. 

nomer in the MIR. Each of our radio-compact objects does indeed 
have a MIR-bright core with apparent scale less than 1'. Some 
are well separated from the surrounding ISM (G31 1.42+0.60, Cs: 
Fig. ll6t . but there are often associated multiple centres of emission, 
intersecting filaments (G31 1.59-0.59, Cf: Fig.HTt. and even ex- 
tended haloes in which 8.0-/^m emission dominates. G310. 89+0.01 
(Cs: Fig. ll8i has a halo in the form of an extensive streamer but its 
most obvious feature in the MIR is the bright partial shell around 
its core. G312.60+0.06, (Ch: Fig.[T9l is a particularly clear case of 
a halo around a compact source. It is associated with a small white 
core (implying MIR thermal continuum) but it is flanked by quite 
extensive and relatively bright diffuse emission. The haloes could 
be emitting in PAH bands, H2 lines or simply be due to contin- 
uum from warm (~350 K) dust grains heated by starlight. The cores 
of compact Hll regions are generally yellow (PAHs) but some are 
white, indicative of the role played by thermal emission from hot 
dust rather than PAHs or any specific fine structure or H2 emission 
lines. 

Most shell-type Hll regions consist of bounded ionized zones 
in the radio and in Ha emission. Their MIR counterparts are fil- 
amentary PDRs that enclose the ionized gas. Some may include 
dominant, bright cores (G31 1.48+0.37, Cs: Fig.l20l 

Our sample contains relatively few shell sources. 
G3 12.68+0.04 (type Si: Fig. |2l} is the north-western of a 



22 Martin Cohen et al. 




G31 1.42+0.60 




Figure 16. Compact Hll region G3 1 1 .42+0.60 (Cs) with shell and multiple 
cores. Image is 162 by 288 arcsec. 




.G31 1.59-0.59 - 



Figure 17. Compact Hll region G31 1.59—0.59 (Cf), with its filamentary 
core. Image is 169 by 240 arcsec. 

pair of interacting bubbles. The second object is G3 12.7 1+0.02 
(Si), south-east of the bright PDR that arises from the interaction 
of the shells blown by the winds from the two exciting stars. The 
odd structure, G312.67— 0.12 (S) (Fig. 1221 consists of a bundle 
of filaments with contiguous diffuse emission, suggestive of a 
partial shell. As noted by CG, even at the resolution of MSX, the 
shell comprising the PDR of G310.99+0.42 (S: Fig.|23j is entirely 
composed of overlapping filaments. These are very well shown at 
the higher resolution of SST (see also Fig.0. 



Figure 18. Compact Hll region G310.89+0.01 (Cs), with an extended halo 
in the form of a streamer, as well as a bright, shell at its core. Image is 240 
by 240 arcsec. 



G312.60+0.J05 




Figure 19. Compact Hll region G312.60+0.05 (Ch), with a blight, diffuse 
halo. Image is 240 by 240 arcsec. 

Diffuse regions are the second most numerous morphology 
in the MIR. Truly diffuse emission (no filaments, shell, or core) 
is associated with G3 11.50-0.48 (type D:Fig.|24]. Most of their 
MIR emission is diffuse and purple, a colour that might indicate 
that H2 lines dominate both the 4.5 and 8.0-/im bands. Its IRAS 
LRS spectrum (Fig.|5j is both too noisy and at too poor a spectral 
resolution to be capable of showing any of these molecular lines. 
Diffuse regions may contain filaments, such as G313.07+0.32 (Df: 
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Figure 20. Compact Hll region G31 1.48+0.37 (Cs) which is dominated by 
a bright compact region yet incorporates both shells and filaments in the 
MIR. Image is 720 by 480 arcsec. 



G312.68+0.04 




Figure 21. Pair of interacting shell regions with bright interaction bar: 
G312.68+0.04 (Ci) and G312.71+0.02 (Ci). Image is 360 by 480 arcsec. 

Fig.|25j. Its MIR emission is red because it is detected only at 8.0- 
/im. Fig. l25l is a fine example of an object with so much diffuse PAH 
emission that it would be almost impossible to determine where the 
real sky background occurs within this box, in order to set it to a 
black level. 

G312.45+0.08 fFig. !26> is a spectacular network of filaments 
over 1° in length, containing a pair of braided filaments as if two 
very large bubbles have come into contact without disruption. The 




Figure 22. Partial shell source. G3 12.67-0. 12 (S). Image is 288 by 
240 arcsec. 
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Figure 23. Shell Hll region G3 10.99+0.42 (S). Image is 480 by 720 arcsec. 
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Figure 24. Diffuse Hll region G31 1.50-0.48 (D). Image is 720 by 
720 arcsec. 
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Figure 25. Diffuse Hll region, G3 13.07+0.32 (Df), containing filaments. 
Image is 480 by 360arcsec. 




Figure 26. Large filamentary region G3 12.45+0.08 (F). Image is 70 by 
60 arcmin. 

interaction is apparently not vigorous but has led to a warping of 
the interface between them (CG). 

Hll regions with compact MIR morphology all have compact 
(or ultracompact) radio morphology. Filamentary MIR regions are 
likewise filamentary radio regions. But shells and diffuse regions 
show no unique relationship between radio and MIR morphologies, 
most likely because the resolution of the SST is so much higher 
than that of the MOST, but also because PAHs arise in the neutral 
zones that surround Hll regions, giving them complementary, rather 
than identical, structures. 

SST's high resolution implies that one could characterize Hll 
regions solely using their MIR structure. Kurtz et al. (1999) demon- 
strated that, in the sample of 15 Hll regions they studied in the radio 
continuum, extended emission is found in 12 objects. The dynamic 



range required to detect these haloes can be assessed from their 
maps given that their contour levels correspond to factors of 2, the 
lowest levels are given, as are the rms noise values in their Very 
Large Array (VLA) maps. Typical dynamic ranges are between 60 
and 500, calculated as the ratio between the compact peak level 
and the 3-4 rms contour. The integrated extended emission in a 
5 arcmin square centred on the compact components ranges from 
1 to 15 times that of the peak intensity, in those objects for which 
Kurtz et al. (1999) regard a physical association between halo and 
compact sources to be possible. One third (six regions) of our MIR 
compact sources show extended emission. The dynamic range in 
surface brightness from each 8.0-/im peak to the faintest extended 
emission clearly detected with 1RAC (all quantities are measured 
above the sky emission) ranges from 8 to 50. The integrated 8.0-/^m 
brightness of these haloes, compared with a point source matched 
to the compact peaks' surface brightness levels, lies between 10 
and 140. Thus, the MIR dynamic range probed by the GLIMPSE 
survey is about 10 times smaller than that at 3.6 cm (the wavelength 
used by Kurtz and his colleagues). If associated with the compact 
regions, the MIR haloes would contribute 10 times more flux rela- 
tive to their compact peaks than at 3.6 cm. If the MIR haloes radi- 
ated thermally then the expected 843-MHz continuum from these 
haloes would be 36 times smaller than that seen by SST (using our 
median ratio of IRAC/MOST without recalibration). The 3.6 cm ra- 
dio emission would be 47 times smaller than at 8.0 /im, with the 
frequency dependence of the Gaunt factor for a typical Hll re- 
gion. The average rms noise in the VLA survey fields (Kurtz et al. 
1999, their Table 2) is 0.5 mjy beam -1 so the 3 rms level is about 
1.5 mjy beam -1 , equivalent to 38 MJy sr -1 at 8.0 [im. The average 
of the 8.0-/im sky-subtracted halo brightnesses for the six regions 
with MIR type "Ch" is 35±10MJy sr~\ Therefore, GLIMPSE at 
8.0 /im provides comparable sensitivity to extended MIR emission 
around compact Hll regions to that afforded by VLA surveys such 
as that by Kurtz and colleagues. 

Are there MIR haloes akin to the extended weakly ionized 
radio haloes? Within the ionized gas of an Hll region PAHs suf- 
fer photodestruction by UV photons. If any survived, the wind of 
the central star would rapidly sweep them away into the PDR. Al- 
lain, Leach & Sedlmayr (1996a) have shown that small PAHs (those 
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with fewer than 50 carbon atoms) are destroyed by photodissocia- 
tion on a timescale of a few years in those places in which the char- 
acteristic MIR bands are seen, specifically for Hll regions. Cations 
and partially hydrogenated PAHs may survive for only about a year 
(Allain, Leach & Sedlmayr 1996b) because they are less stable than 
the parent molecules. Consequently, small PAHs must be made in 
abundance in locations in which the MIR bands are emitted, for 
example in shocks with speeds between 50 and 200 kms -1 where 
shattering and fragmentation create PAH molecules, PAH clusters, 
and small grains (Jones, Tielens & Hollenbach 1996). PAHs with 
more than 50 carbon atoms may survive the UV radiation field 
and even small PAHs might last as long as 1000 years in the dif- 
fuse ISM. Therefore, PAHs might exist for significant periods in 
extended haloes that were only weakly ionized and far from the 
dominant source of destructive UV photons. If all the MIR halo 
emission were due to PAHs then it would be plausible to associate 
weak diffuse PAH emission with extended radio haloes in which 
low fluxes of UV photons leaking from dense ionized cores excite 
rather than destroy PAHs. MIR spectroscopy of such haloes would 
be definitive. 



9 G311.70+0.31: A UNIQUE STRUCTURE 

Fie. 1271 presents the point-source-subtracted image of a most un- 
usual structure in this field, weakly detected by MSX, but clearly 
present in all four IRAC bands. Its spine consists of a series of 
segments that strongly suggests that a Rayleigh-Taylor instability 
might be responsible. If so then the mean wavelength of the seg- 
ments is 38±3" (standard error of the mean). However, lacking any 
estimate of distance and optical forbidden-line spectroscopy, we are 
unable to use this information on what is surely the fastest growing 
Rayleigh-Taylor mode to derive a consistent set of physical param- 
eters for the nebula. 843-MHz radio emission was detected by the 
MOST and is displayed by white contours over the SST 8.0-/im im- 
age. The object's ratio of MSX 8.3-^im integrated flux to that from 
the MOST at 843 MHz is 27±2, entirely consistent with thermal 
radio emission and PDRs associated with PAHs. The presence of 
PAHs is confirmed by the identical morphologies in the IRAC 5.8 
and 8.0-/im bands. The radio emission has two peaks, one within 
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Figure 27. Upper: G31 1.70+0.31 IRAC 8.0-/im residual image. Circular 
artefacts are the remnants of point sources that have been subtracted from 
the original image, in particular a very bright object at (311.704,0.294). 
Lower: G31 1.70+0.31 8.0-/im image overlaid by 843-MHz contours at lev- 
els of 2, 5, 7, 10, 12 mJy beam - 1 . 



the "bay" of MIR emission to the east, and a second comparable 
peak west of the centre. A physical association with the east peak 
is morphologically stronger. If the source represented a structure 
blown and ionized by a stellar wind then one would expect the hot 
star to betray its presence by a radio continuum peak and by MIR or 
far-infrared emission. A search of the region indicates only a single 
plausible external source that might be a star with a wind, namely 
/AA514028— 6103, a 25-/im-only object. Nothing definitive can be 
gleaned from the IRAS photometry or its Low Resolution Spec- 
trum. No counterparts appear in either the 2MASS images or those 
of the Southern Ha Survey (Parker et al. 2005), suggesting high 
extinction toward the putative MIR and radio source. 



10 GALAXIES, SNRS, AND OTHER NONTHERMAL 
EMITTERS 

GLIMPSE at 8.0 /im offers a survey three orders of magnitude 
deeper than MSX at 8.3 /im. This enables us to probe the MIR/radio 
ratios for sources other than Hll regions. All five of the sus- 
pected background galaxies in the /=312° field, described by CG, 
were detected by the SST at 8.0 fim. All appear as point sources 
and so require no aperture correction. The ratios of IRAC/MOST 
for G3 10.86+0.01, G3 1 1 .54-0.02, G3 12. 1 1 -0.20, G3 13.25+0.32, 
and G313.42+0.09 are 0.080, 0.009, 0.081, 0.019, and 0.034, re- 




Figure 29. SNR G311. 52-0.34 false colour image: blue 3.6, green 4.5, red 
5.8 /im. The SNR corresponds to the complete ring around the centre. 



spectively. The average ratio for these galaxies is 0.04±0.01, 900 
times smaller than the empirical ratio for thermal emission. 

If we interpret the extended 8-/im structure that is associated 
with radio thermal emission as due to PAHs then it is interesting 
that even diffuse radio regions in this field have 8-pim emission. 
This implies that these patchy zones are not optically thin to ioniz- 
ing radiation because they apparently have PDRs. Seemingly the 
only environment in this field that shows a deficit of PAHs (as 
traced by 8-/im imagery) is toward SNRs. 

Very few SNRs were detected by MSX (Cohen et al. (2005) 
present an MSX image of Cas-A) and nonthermal emission regions 



are very faint below 9 /J,m, compared with thermal zones. When 
one might suspect that a nonthermal region was detected it is in- 
variably the case that MSX observed a mixture of thermal and non- 
thermal filaments but only thermal structures were detected. Two 
(G310.65-0.29, G310.80-0.41) of the four SNRs in this field are 
so confused by diffuse and filamentary thermal emission that no 
useful measurements or even upper limits can be made of their 
purely nonthermal portions. The SNR G3 12.44—0.36 is very large 
(40 arcmin) and its interior includes a considerable amount of dif- 
fuse emission, probably thermal, and a considerable number of 
point and compact sources. These circumstances render it a poor 
candidate for determining the ratio of nonthermal MIR/radio emis- 
sion. 

However, the radio-bright young remnant G31 1.52— 0.34, dis- 
covered by Green (1974), is detected by GLIMPSE. Fig.l28lover- 
lays the 4.5-/im residual IRAC image by 13-cm continuum contours 
(ATCA data kindly provided by Josi Gelfand). The spatially inte- 
grated, sky-subtracted 8.0-/im flux density of this SNR can be com- 
pared with the integrated radio continuum. The empirical 8.0-/^m 
to 843-MHz ratio is 0.19, or 0.14 after aperture correction. WCL 
list one other nonthermal object, G313.45— 0.13, but this too lies 
amidst bright thermal filaments. Combining the ratios for the five 
galaxies with the aperture-corrected ratio for SNR G31 1.52—0.34 
yields an average IRAC/MOST ratio of 0.06±0.02, over 600 times 
smaller than the thermal ratio. On this basis, the MIR/radio ratio is 
a sharp discriminant between thermal and nonthermal emission. 

The electronic version of our paper carries a 3-band false 
colour image (Fig. I29> in which 3.6, 4.5, 5.8-/im emission are 
coded blue, green and red, respectively and the remnant can be 
traced around its complete shell. We have sampled the four bright- 
est areas of MIR emission (roughly one in each quandrant around 
the elliptical ring) to derive four estimates of IRAC colours for the 
SNR. The mean values that we calculate yield observed ratios of 
integrated flux densities at 3.6 and 5.8 /im of 0.17±0.06, and at 4.5 
and 8.0 fim of 0.43±0.06. After applying the SSC aperture correc- 
tions these become 0.21±0.07 and 0.55±0.08, respectively. Reach 
et al. (2006) suggest that shocked gas contributes the IRAC emis- 
sion from this SNR. Our colours, within their uncertainties, overlap 
the predicted colour-colour box presented by Reach et al. (2006: 
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their Fig. 2) for ionized gas that has suffered a fast shock. The MIR 
emission is most obvious at 4.5 /im, and somewhat less bright at 
5.8 pim, due dominantly to Br-a and [Fell] lines, respectively, ac- 
cording to Fig. 1 of Reach et al. (2006). 



11 CONCLUSIONS 

We present evidence that the absolute diffuse calibration of SST 
at 8.0 /im is too high by about 36 percent. All measured IRAC sur- 
face brightnesses of any diffuse emission region should be scaled 
down by a multiplicative factor of 0.74±0.07, based on the abso- 
lutely calibrated MSX images. This provides a totally independent 
confirmation of the separate analyses of IRAC's diffuse calibration 
that are based upon the fundamental calibration of the instrument 
(0.737: Reach et al. 2005) and detailed modeling of the light distri- 
butions of elliptical galaxies (0.74±0.07: Jarrett 2006). Therefore, 
the recommendation of the SSC that one should scale down 8.0-/im 
diffuse BCD flux densities by a multiplicative factor of 0.74 should 
indeed produce the best absolute estimate of any extended 8.0-/im 
emission. 

We find that MIR and radio fluxes are linearly related from 
our regression analysis. These correlations reflect the intimate as- 
sociation between ionized zones and their surrounding PDRs in 
which polycyclic aromatic hydrocarbon band emission arises. Our 
median analysis for Hll regions of Fg.3 Mm /Ss43 MHz yields 25±5 
(standard error of the median), confirming CG's original deriva- 
tion. This ratio of MIR/radio integrated fluxes is independent of 
morphology of the ionized region but is valid only for radio struc- 
tures smaller than ~24 arcmin. Aperture synthesis by MOST im- 
plies that the telescope is progressively less sensitive to structure on 
scale sizes >15 arcmin, diminishing the estimated fluxes of larger 
sources. Compact Hll regions below ~3 arcmin in size tend to have 
somewhat smaller MIR/radio ratios, which might be attributed to 
absorption by dust grains of UV photons in some of these small, 
dense regions. The only cometary region in our sample appears to 
suffer substantial thermal self-absorption of its radio flux. 

We find an equivalent correlation for Hll regions between ob- 
served 8.0 /im and 843-MHz fluxes, albeit with more scatter than 
for MSX. With the SST pipeline and calibration that were applied 



to the GLIMPSE images released to the community, we find an 
observed value of F8.0Mm/S843 mhz = 36±13 best represents this 
relationship. Removing the 36 percent absolute calibration problem 
this corresponds to an aperture-corrected ratio for IRAC 8.0/MOST 
of 27±10, very similar to the MSX/MOST criterion for radio ther- 
mal emission regions. Supernova remnants studied with GLIMPSE 
images (Reach et al. 2006) show that none has appreciable MIR 
synchrotron radiation and the dominant post-shock cooling is usu- 
ally through H2 lines. This leads to much smaller MIR intensities 
compared with the radio continuum than for thermal sources, which 
are intimately associated with PDRs and their bright and broadband 
PAH emission features. The MIR/radio flux ratio discriminates be- 
tween thermal and nonthermal regions. For the latter, the ratio is 
0.06±0.02, some 600 times smaller than for Hll regions. 

The MIR morphology of Hll regions does not match their ra- 
dio structure one-to-one. This is most readily seen by examina- 
tion of Table 1 where the new homogeneous radio and MIR mor- 
phologies appear side by side. A diversity of MIR characteristics 
can be associated with individual radio types. There is clear evi- 
dence for extended MIR emission surrounding even compact and 
cometary Hll regions. Either filaments or diffuse emission may be 
found around the cores of compact regions, though some appear 
truly isolated on the sub-arcminute scale. When diffuse MIR struc- 
ture is found outside the cores of compact Hll regions it has the 
colour of PAH emission. Consequently, compact Hll regions may 
indeed have MIR counterparts to the postulated weakly ionized ra- 
dio haloes. A valid question after examination of these GLIMPSE 
images is whether there are any compact Hll regions in the MIR 
without shells, filaments, and haloes. Given the high-resolution of 
the SST with IRAC one might wish to develop new quantitative 
criteria for Hll regions, in addition to our MIR morphological char- 
acterization. 

The SEDs of all the Hll regions in our sample show little vari- 
ation in the wavelength range seen by the four IRAC bands. Well- 
defined colours have been determined from this sample, and diag- 
nostic boxes proposed for two IRAC colour-colour planes. There is 
good separation in these plots between the colours of diffuse ther- 
mal emission and those of the plethora of MIR point sources en- 
countered in the Galactic Plane. The colours of known pointlike or 
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ultracompact Hll regions are also sharply distinguished from those 
of diffuse thermal regions in these plots. 

The advantages of combining MIR and radio aspects of par- 
ticular sources are illustrated by the unique MIR structure of 
G3 1 1.70+0.3 1 in this field. IRACs two longest bands offer a highly 
sensitive tool for testing for PAH emission by comparing mor- 
phologies at different wavelengths that are each dominated by these 
molecular bands. False colour imagery of G31 1.70+0.31 shows the 
excellent match between the MIR configurations of the source at 
5.8 and 8.0 /im, strongly suggestive of the dominance of PAH spec- 
tral features in these two images. The flux ratio, MSX/MOST, is ex- 
actly what one expects for the combination of thermal radio contin- 
uum and 8.3-/im emission from PAH features. CG emphasized the 
value of MSX's filters in which PAH structure detected at 8.3 /J,m 
would be mimicked in the 12.1-/im image. MSXhad a large dispar- 
ity in responsivity between the 8.3-^im band and all others, whereas 
IRAC offers four sensitive detectors so false colour imagery can 
achieve its true potential. 
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